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Cranked Relativistic Hartree-Bogoliubov (CRHB) theory is presented as an extension of Relativistic 
Mean Field theory with pairing correlations to the rotating frame. Pairing correlations are taken 
into account by a finite range two-body force of Gogny type and approximate particle number 
projection is performed by Lipkin-Nogami method. This theory is applied to the description of 
yrast superdeformed rotational bands observed in even-even nuclei of the A ~ 190 mass region. 
Using the well established parameter sets NLl for the Lagrangian and DIS for the pairing force one 
obtains a very successful description of data such as kinematic (J^'^') and dynamic (J''^-') moments 
of inertia without any adjustment of new parameters. Within the present experimental accuracy 
the calculated transition quadrupole moments Qt agree reasonably well with the observed data. 

PACS numbers: 21.60.-n, 21.60.Cs, 21.60. Jx, 27.80.-Fw 



The investigation of superdeformation in different mass 
regions still remains in the focus of low-energy nuclear 
physics. Experimental data on superdeformed rotational 
(SD) bands are now available in different parts of the pe- 
riodic table, namely, in the >1 ~ 60 80, 130, 150 and 
190 H] mass regions. This richness of data provides the 
necessary input for a test of different theoretical models 
and the underlying effective interactions at superdefor- 
mation. Cranked relativistic mean field (CRMF) theory 
developed in Refs. [||-|5| represents one of such theories. 
It has been applied in a systematic way for the descrip- 
tion of SD bands observed in the A ^ and A ~ 150 
mass regions. The pairing correlations in these bands are 
considerably quenched and at high rotational frequencies 
a very good description of experimental data is obtained 
in the unpaired formalism in most of the cases as shown 
in Refs. 

On the contrary, pairing correlations have a consider- 
able impact on the properties of SD bands observed in the 
A ~ 190 mass region and more generally on rotational 
bands at low spin. Different theoretical mean field meth- 
ods have been applied for the study of SD bands in this 
mass region. These are the cranked Nilsson-Strutinsky 
approach based on a Woods-Saxon potential |10 1^, self- 
consistent cranked Hartree-Fock-Bogoliubov approaches 
based either on Skyrme |l^,|l^ or Gogny force s ||l^,[l5|| . It 
was shown in different theoretical models (ll - 14 16|] that 
in order to describe the experimental data on moments of 
inertia one should go beyond the mean field approxima- 
tion and deal with fluctuations in the pairing correlations 
using particle number projection. This is typically done 
in an approximate way by the Lipkin-Nogami method 
-M. With exception of approaches based on Gogny 



forces, special care should also be taken to the form of the 
pairing interaction. For example, quadrupole pairing has 
been used in addition to monopole pairing in the cranked 
Nilsson-Strutinsky approach jl|]. A similar approach to 
pairing has also been used in projected shell model [po[ . 
Density dependent pairing has been used in connection 
to Skyrme forces |12|. These requires, however, the ad- 
justment of additional parameters to the experimental 
data. 

Cranked Relativistic Hartree-Bogoliubov (CRHB) the- 
ory presented in this article is an extension of cranked 
relativistic mean field (CRMF) theory to the description 
of pairing correlations in rotating nuclei. A brief outline 
of this theory and its application to the study of sev- 
eral yrast SD bands observed in even-even nuclei of the 
A ~ 190 region with neutron numbers N = 110, 112, 114 
is presented below while more details (both of the the- 
ory and the calculations) will be given in a forthcoming 
publication. 

The theory describes the nucleus as a system of Dirac 
nucleons which interact in a relativistic covariant manner 
through the exchange of virtual mesons : the isoscalar 
scalar a meson, the isoscalar vector uj meson, and the 
isovector vector p meson. The photon field (A) accounts 
for the electromagnetic interaction. 

The CRHB equations for the fermions in the rotating 
frame are given in one-dimensional cranking approxima- 
tion by 
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where ft. = /id — A is the single-nucleon Dirac Hamilto- 
nian minus the chemical potential A and A is the pairing 
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potential. Jx and 0,^ are the projection of total angular 
momentum on the rotation axis and the rotational fre- 
quency. Uk and Vfc are quasiparticle Dirac spinors and Ek 
denote the quasiparticle energies. The variational prin- 
ciple leads to time-independent inhomogeneous Klein- 
Gordon equations for the mesonic fields in the rotating 
frame 



the spatial components of the vector potential A{r) are 
neglected. 
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where the source terms are sums of bilinear products of 
baryon amplitudes 
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The sums over fc > run over all quasiparticle states cor- 
responding to positive energy single-particle states {no- 
sea approximation). In Eqs. (^,^, the indexes n and 
p indicate neutron and proton states, respectively, and 
the indexes is and iv arc used for isoscalar and isovector 
quantities. Pyir), f\r) in Eq. (||) correspond to p^^ir) 
and j^'^{r) defined in Eq. (||), respectively, but with the 
sums over neutron states neglected. 

The spatial components of the vector mesons give ori- 
gin to a magnetic potential V(r) which breaks time- 
reversal symmetry and removes the degeneracy between 
nucleonic states related via this symmetry . This ef- 
fect is commonly referred as a nuclear magnetism . It is 
very important for a proper description of the moments 
of inertia [Q. Consequently, the spatial components of 
the vector mesons w and p are properly taken into ac- 
count in a fully self-consistent way. Since the coupling 
constant of the electromagnetic interaction is small com- 
pared with the coupling constants of the meson fields, 
the Coriolis term for the Coulomb potential AQ{r) and 
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FIG. 1. Dynamic (J'^') (upper panels) and kinematic 
(j'^^) (bottom panels) moments of inertia. The experimen- 
tal values (full circles) are compared with CRHB calculations 
(lines). Solid and dashed lines show the results of calcula- 
tions with and without approximate particle number projec- 
tion performed by means of Lipkin-Nogami method, respec- 
tively. The experimental data are taken from Refs. |26|-|2^] 
(i«4pb) and (is^Hg). 

In the present version of CRHB theory, pairing corre- 
lations are only considered between the baryons, because 
pairing is a genuine non-relativistic effect, which plays a 
role only in the vicinity of the Fermi surface. The phe- 
nomenological Gogny-type finite range interaction 

,-[(ri-r2)M]' 
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with the parameters pi, Wi, Bi, Hi and Mi (i — 1,2) 
is employed in the pp (pairing) channel. The parameter 
set DIS [|2| has been used in the present calculations. 
This procedure requires no cutoff and provides a very 
reliable description of pairing properties in finite nuclei. 
In conjuction with relativistic mean field theory such an 
approach to the description of pairing correlations has 
been applied, for example, in the study of ground state 
propertie s [p5[ , neutron halos p^ , and deformed proton 
emitters [p4[. In the present approach we go beyond the 
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mean field and perform an approximate particle num- 
ber projection before the variation by means of Lipkin- 
Noeami method As illustrated in Fig. 0, this 

feature is extremely important for a proper description 
of the moments of inertia. 
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FIG. 2. The same as in Fig. |l| but for the yrast SD bands 
in I92,i96p|_| rpj^g experimental data are taken from Refs. ^] 
(^^^Pb) and [@ ("'^Pb). The experimental j'^' values of the 
^^^Pb band are shown for the spin value of Jo = 8 for the 
lowest state in the SD band (the lowest transition with the 
energy 214.8 keV corresponds to a spin change 10 ^8). The 
experimental J^^^ values of the ^^^Pb band are shown for three 
different values of /q. The values being in best agreement with 
calculations are indicated by solid circles. 

The present calculations have been performed with the 
NLl parametrization |^ of the relativistic mean field 
Lagrangian. The CRHB-equations are solved in the basis 
of an anisotropic three-dimensional harmonic oscillator in 
Cartesian coordinates. A basis deformation of f3o = 0.5 
has been used. All fermionic and bosonic states belonging 
to the shells up to A'^ = 14 and Nb — 16 are taken into 
account in the diagonalisation and the matrix inversion, 
respectively. This truncation scheme provides reasonable 
numerical accuracy. For example, the increase of the 
fermionic basis up to Np — 17 changes the values of 
the kinematic moment of inertia j'^' and the transition 
quadrupole moment Qt by less than 1%. The numerical 
errors for the total energy are even smaller. 
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FIG. 3. The same as in Fig. ^ but for the yrast SD bands 
in ^^°'^^^Hg. The experimental data are taken from Refs. ^] 
(i«°Hg) and @ ("^Hg). 

Yrast SD bands in ^^''Pb and ^^"'Hg are linked to the 
low-spin level scheme |2^,^,Q. In addition, there is a 
tentative hnking of SD band in ^^^Pb Q. These data 
provide an opportunity to compare with experiment in 
a direct way not only calculated dynamic (J^^^) but also 
kinematic (J^^^) moments of inertia. On the contrary, at 
present the yrast SD bands in i90,i92jjg ^^^^ ^^^Ph are 
not linked to the low-spin level scheme yet. Thus some 
spin values consistent with the signature of the calculated 
yrast SD configuration should be assumed for the exper- 
imental bands when a comparison is made with respect 
of the kinematic moment of inertia J^^\ 

The results of such a comparison are shown in Figs. 
|l|, I and |. The theoretical values agree well with 
the experimental ones in the cases of linked SD bands 
in ^^''Pb and ^^"^Hg and tentatively linked SD band in 
^^^Pb. The comparison of theoretical and experimental 
J*^^^ values (see Figs. || and ||) indicates that the low- 
est transitions in the yrast SD bands of ^^°Hg, ^^^Hg 
and i^spb with energies 316.9, 214.4 and 171.5 keV, re- 
spectively, most likely correspond to the spin changes of 
14+ ^ 12+, 10+ 8+ and 8+^6+. If these spin values 
are assumed, good agreement between theory and exper- 
iment is observed. Calculated and experimental values of 
the dynamic moment of inertia J^^^ agree also well, see 
Figs. |1|, I and |. 

The increase of kinematic and dynamic moments of in- 



3 



ertia in this mass region can be understood in the frame- 
work of CRHB theory as emerging predominantly from 
a combination of three effects: the gradual alignment 
of a pair of J15/2 neutrons, the alignment of a pair of 
*i3/2 protons at a somewhat higher frequency, and de- 
creasing pairing correlations with increasing rotational 
frequency. The interplay of alignments of neutron and 
proton pairs is more clearly seen in Pb isotopes where 
the calculated J^^^ values show either a small peak (for 
example, at ~ 0.45 MeV in ^^^Pb, see Fig. |) or a 
plateau (at ~ 0.4 MeV in ^^^pb^ gee Fig. |). With in- 
creasing rotational frequency, the J^^-* values determined 
by the alignment in the neutron subsystem decrease but 
this process is compensated by the increase of due 
to the alignment of the 113/2 proton pair. This leads to 
the increase of the total j(2)-value at > 0.45 MeV. 
The shape of the peak (plateau) in J*^^^ is determined 
by a delicate balance between alignments in the proton 
and neutron subsystems which depends on deformation, 
rotational frequency and Fermi energy. For example, 
no increase in the total dynamic moment of inertia j'^' 
has been found in the calculations after the peak up to 
fl^ = 0.5 MeV in ^^^Hg, see Fig. ^. It is also of interest 
to mention that the sharp increase in j'^' of the yrast 
SD band in ^^^Hg is also reproduced in the present cal- 
culations. One should note that the calculations slightly 
overestimate the magnitude of J*^^) at the highest ob- 
served frequencies. The possible reasons could be the 
deficiencies either of the Lipkin-Nogami method jssj or 
the cranking model in the band crossing region or both 
of them. 

The comparison between calculated and experimen- 
tal absolute transition quadrupole moments Qt is less 
straightforward. This is because the uncertainties in ab- 
solute measured Qt values arising from the uncertainties 
in stopping powers can be as large as 15% Thus 
the comparison of Qt's values obtained in different ex- 
periments should be performed with some caution since 
systematic errors due to different stopping powers may 
be responsible for the observed differences. In addition, 
as illustrated in Fig. ^, the experimental Qt values de- 
pend somewhat on the type of analysis (centroid shift or 
line shape) used when these quantities are extracted from 
the data. 

The results of CRHB calculations are compared with 
most recent experimental data in Fig. ^. One can con- 
clude that the results of calculations for absolute values 
of Qt will be within the 'full' error bars if the 15% un- 
certainty due to stopping powers would be taken into 
account (experimental data shown in Fig. ^ does not in- 
clude these uncertainties). For the sake of simplicity we 
will not take into account these uncertainties in the sub- 
sequent discussion and will concentrate mainly on the 
experimental data obtained with the same stopping pow- 
ers. In Fig. H such data are indicated by the same cap- 



ital letters. While the calculated Qt values are close to 
the experimental values obtained with centroid shift and 
line shape analysis for Pb isotopes, most of experimen- 
tal Qt values (with exception of exp. B) are overesti- 
mated in calculations in the case of Hg isotopes. One 
should note that the most recent experimental data on 
i92jjg jg contradictory since two experiments (exp. A 
p7t and exp. B iQ) give very different values of Qt, see 
Fig. 1^. Definitely, the measurements of relative transi- 
tion quadrupole moments between SD bands in Pb and 
Hg isotopes using the same stopping powers, which are 
not available nowadays, are needed to find out, whether 
this discrepancy between calculations and experiment is 
due to an inadequate theoretical description or the ex- 
perimental problems quoted above. In the calculations, 
relative average quadrupole moments AQt between yrast 
SD bands of Pb and Hg isotopes decrease with increasing 
neutron number N {AQt ~ 1.6 eb, « 1.4 eb and « 1.06 
eb for N = 110, 112 and 114, see Fig. |). 
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FIG. 4. Measured and calculated transition quadrupole 
moments Qt of the yrast SD bands. Solid and dashed lines 
show the results of calculations for Hg and Pb isotopes, re- 
spectively. Experimental values are given by error bars where 
the type of symbols corresponding to experimental middle 
value of Qt indicates the isotope and type of analysis used to 
extract the experimental value. The following correspondence 
exists: circles - Pb isotopes, triangles - Hg isotopes, open sym- 
bols - line shape analysis, filled symbols - centroid shift. The 
experimental data are taken from Refs. (exp. A), ^] 
(exp. B), (exp. C), and (exp. D). Experimental data 
on 
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Results of calculations indicate the general trend of de- 
crease of average Qt values with the increase of neutron 
number N both for Pb and Hg isotopes. The results of 
the centroid shift analysis for i94,i96p|-, (^g^p. C) indi- 
cate a slight decrease in the Qt values with increasing N 
consistent with theoretical results. Although the data on 
I92,i94jjg (^gxp. A) indicate similar values of Qt being in 
slight contradiction with theoretical results, definite con- 
clusions are not possible at present due to the large error 
bars. In addition, with increasing rotational frequency 
U,x the calculated Qt values show initially a slight in- 
crease which is followed by a subsequent decrease. In the 
case of i^^Hg this feature is hidden by the band crossing. 
The maximum Qt values within specific configurations 
are calculated at different frequencies Vlx as a function of 
the neutron number TV. With increasing N the maximum 
Qt is reached at lower frequencies. Similar variations of 
Qt have been also observed in the cranked Hartree-Fock 
calculations with Skyrme forces Dedicated experi- 

ments aiming on the measurements of the variations of 
transition quadrupole moments Qt as a function of ro- 
tational frequency are needed in order to confirm or 
reject these results. 

In conclusion, the Cranked Relativistic Hartree- 
Bogoliubov theory has been developed and applied to the 
description of yrast SD bands observed in the A ~ 190 
mass region. With an approximate particle number pro- 
jection performed by the Lipkin-Nogami method, the 
rotational features of experimental bands such as kine- 
matic and dynamic moments of inertia are very well de- 
scribed in the calculations. Calculated values of transi- 
tion quadrupole moments Qt are close to the measured 
ones, however, more accurate and consistent experimen- 
tal data on Qt is needed in order to make detailed com- 
parisons between experiment and theory. 
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